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ABSTRACT 

The  sputtering  of  Cu  fee  (110),  (100)  and  (111)  surfaces  by  1- 
to  10-keV  A  ions  at  normal  incidence  has  been  investigated  by  high- 
speed digital  computer  techniques.  The  interatomic  repulsions  are 
described  by  Born-Mayer  potentials.  The  shape  of  the  sputtering 
ratio  curves  are  in  close  agreement  with  the  data  of  Magnuson  and 
Carls ton.  Ejection  patterns  in  accord  with  experimental  observations 
are  predicted.  There  is  evidence  that  the  sputtering  mechanism  is 
not  the  result  of  momentum  focusing. 
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1 .  Introduction . 

The  process  of  ionic  erosion  of  a  solid  surface  is  known  as 
sputtering.  As  pointed  out  by  Cheney,  Rogers  and  Pitkin  fl]  and 
Harrison  (_2J  >  various  theoretical  models  of  the  sputtering  process 
have  been  developed  since  the  1850ts  when  metal  from  the  electrodes 
in  gas  discharge  tubes  was  first  observed  to  be  transferred  to  the 
glass  walls  of  the  tubes.  In  1923,  Kingdon  and  Langmuir  (Vj 
suggested  a  momentum  transfer  model,  later  elaborated  upon  by  Lamar 
and  Corapton  QTJ  .  This  theory  was  subsequently  replaced  by  a 
thermal  spike  model,  wherein  a  local  hot  spot  caused  evaporation 
(5,6J.  Later  Keywell  [7J  suggested  a  model  based  upon  a  succession 
of  random  collisions,  while  Harrison  |_8-13 J  has  proposed  several 
models,  including  a  statistical  approach  based  on  neutron  diffusion 
theory,  and  Monte  Carlo  techniques.  Pease  [14]  developed  a  model 
based  upon  expressions  from  radiation  damage  theory. 

However,  these  early  models  neglected  the  ordered  arrangement 
of  atoms  within  the  crystal  lattice.  Current  theories  are  based 
upon  a  transfer  of  energy  through  collision  processes,  as  developed 
by  Silsbee  fl5 1  and  Liebfried  [l6J.  In  addition  to  conducting  a 
variety  of  experiments,  Harrison  jj.7,18]  combined  these  energy  chain 
concepts  with  previous  statistical  theories  into  a  radiation  damage- 
statistical  model  [2,19j  . 

Recently,  workers  in  the  field  of  radiation  damage  have  employed 
digital  computers  to  make  theoretical  studies.  Robinson  and  Oen  |20j 


have  assumed  that  the  "knock-on"  atom  interacts  through  a  series  of 
binary  (two-body)  collisions,  while  Gibson,  Goland,  Milgram  and 
Vineyard  [2ll ,  Erginsoy,  Vineyard  and  Englert  [_22] ,  and  Gay  and 
Harrison  £23j  have  used  iteration  techniques  to  numerically  integrate 
the  classical  equations  of  motion  of  a  many-body  set  of  atoms. 

The  present  work  uses  the  latter  technique  to  investigate  the 
sputtering  of  three  orientations  of  a  copper  crystal  subjected  to 
bombardment  by  an  argon  ion.  This  particular  combination  of  elements 
was  chosen  for  two  reasons j  first,  there  is  experimental  data  avail- 
able for  comparison  124,25]  and  secondly,  parameters  to  be  used  in 
the  Born-Mayer  potentials  have  been  calculated  (21,26,27] . 

Distinction  should  be  made  between  "sputtering"  and  "transmission 
sputtering."  As  used  in  this  paper,  the  former  pertains  to  atoms 
ejected  from  the  bombarded,  or  front,  surface  of  a  solid,  while  the 
latter  refers  to  atoms  leaving  the  back  surface  of  a  thin  crystal. 
Although  the  possibility  exists  that  an  ejected  atom  will  have  one 
or  more  electrons  stripped  away  in  the  sputtering  process,  this 
study  makes  no  distinction  between  neutral  atoms  and  ions. 

The  sputtering  ratio  is  defined  as  the  number  of  atoms  ejected 
per  incident  ion.  Several  aspects  of  sputtering  have  been  studied 
experimentally  [24, 25,28,29,35],  including  the  sputtering  ratio  as 
a  function  of  primary  (incoming  ion)  energy,  incidence  angle,  crystal 
orientation  and  target  temperature j  the  angular  distribution  of  ejected 
atoms,  and  the  energies  of  these  atoms.  This  study  is  concerned 


quantitatively  only  with  the  sputtering  ratio  as  a  function  of 
mono crystalline  orientation  for  primary  energies  in  the  1-  to  10-keV 
range  and  at  normal  incidence,  and  qualitatively  with  a  description 
of  the  process  whereby  the  crystal  atoms  are  ejected  from  the  lattice, 
and  their  angular  distribution. 


2.  The  Model 

The  model  employed  in  this  computer  simulation  program  is  the 
result  of  the  combination  of  many  ideas  from  previous  work  in  the 
field  of  radiation  damage.  An  early  attempt  by  Oen,  Holmes  and 
Robinson  J31|  assumed  binary  collisions  but  excluded  crystal  lattice 
effects.  This  work  suggested  that  neither  the  exponentially  screened 
Coulomb  potential  nor  the  Bohr  hard  sphere  potential  were  sufficiently 
accurate.  A  later  model  used  by  Robinson  and  Oen  f20J  placed  copper 
atoms  on  their  proper  sites  in  a  face  centered  cubic  lattice. 
However,  this  model  also  considered  collisions  as  binary  and  assumed 
that  the  target  atoms  did  not  interact  amongst  themselves.  The  Bohr 
potential  and  Born-Mayer  potential  with  Gibson  2  parameters  £21] 
were  compared.  Gibson  et,:  al.  f^ll  also  studied  effects  in  the  fee 
copper  crystal,  but  considered  collisions  as  n-body  while  comparing 
several  sets  of  the  two  parameters  of  the  Born-Mayer  potential. 
Similar  work  for  bec  iron  was  done  by  Erginsoy  et,  al.  |22j •  More 
recently,  Gay  and  Harrison  Qj3J  employed  a  somewhat  different  n-body 
model.  This  project  uses  the  Gay-Harrison  model  but  with  further 
modifications,  vizs 

(1)  The  stationary  immoveable  atoms  surrounding  the  crystallite 
have  been  removed.  The  object  of  this  study  is  the  ejection  of  atoms 
from  the  surface,  and  also,  the  lattice  is  stable  before  bombardment 
without  these  stationary  atoms. 

(2)  The  crystallite  size  has  been  increased  to  contain  the 


sputtering  event.  The  (110)  sample  contains  115  atoms,  the  (100) 
has  158,  and  the  (111)  contains  150.  The  (110)  sample  corresponds 
to  ten  planes  in  the  x-direction,  five  planes  in  the  y-direction, 
and  nine  planes  in  the  z-direction.  The  (100)  sample  is  7  x  5  x  9, 
while  the  (111)  is  9  x  5  x  17.  The  total  external  dimensions  of 
these  three  crystallites  are  approximately  the  same  (see  appendix  I). 
As  work  progressed  it  was  readily  determined  which  atoms  were 
sputtered  from  the  crystallite.  Thus,  by  judicious  choice  of  impact 
points,  the  crystallite  may  be  decreased  in  size  to  reduce  computer 
running  time.  This  will  be  done  in  future  studies. 

(3)  A  running  report  of  the  motion  of  every  atom  versus  a 
final  status  report  has  been  incorporated. 

(4)  The  impact  area  for  the  (111)  sample  has  been  chosen  to  be 
representative  of  a  "repeat"  crystallographic  distance  rather  than 
just  of  the  first  two  planes.  The  repeat  distance  is  the  number  of 
layers  that  must  be  traversed  to  repeat  the  crystal  configuration 
[24].  Thus  the  projections  of  the  first  two  layers  of  the  (110) 
sample,  first  layer  of  the  (100)  sample,  and  first  three  layers  of 
the  (111)  sample  are  considered.  Then  symmetry  shows  that  the 
rectangle  for  the  (110),  and  triangles  for  the  (100)  and  (111) 
orientations,  will  completely  describe  the  crystallites  if  repeated 
(see  figs.  1-3).  Thirty-six  impact  points  were  assigned  to  the  (110) 
representative  rectangle,  while  21  points  were  assigned  to  the  (100) 
and  (111)  triangles  (see  fig.  4). 


(5)  There  is  no  "target  atom"  as  defined  by  Gay  and  Harrison, 
but  the  impact  points  are  all  near  a  "central"  atom  of  the  crystal- 
lite front  surface. 

(6)  The  "time  step  multiplier"  is  0,10  for  all  impact  points 
except  those  in  which  a  "head-on"  collision  with  a  front  surface 
lattice  atom  occurs.  In  the  latter  case,  a  value  of  0.05  was  found 
to  be  necessary  to  minimize  total  energy  error. 

In  addition  to  n-body  collisions,  major  ideas  used  by  Gay  and 
Harrison  and  retained  in  this  program  include  the  double  iteration 
technique,  lack  of  thermalizing  motion,  concepts  of  representative 
areas  and  time  step  intervals  mentioned  previously,  the  "lattice 
unit,"  disregard  of  interactions  between  the  primary  and  lattice 
electrons  (although  provision  is  made  for  inclusion  in  future  studies), 
use  of  Gibson  2  parameters  in  the  Born-Mayer  potential  for  copper- 
copper  interactions,  and  omission  of  cohesive  effects  and  restoring 
forces.  The  last  point  is  justified  by  time  considerations;  the 
sputtering  event  is  complete  in  a  relatively  short  period,  one-tenth 
to  one-hundreth  of  the  time  necessary  for  these  effects  to  commence 
[23] . 

Harrison,  Carlston  and  Magnuson  |_26j  state  that  the  repulsive 
exponential  Born-Mayer  potential  function  of  the  form  V(r)  =  exp(A+Br), 
where  B  is  always  negative,  appears  to  be  adequate  for  most  high 
energy  collision  dynamics  in  a  lattice.  Initially,  the  program  was 
run  using  the  Gibson  2  parameters  in  the  Born-Mayer  potential  for 


both  the  argon-copper  and  copper-copper  interactions.  These  para- 
meters are  not  valid  for  the  argon-copper  system.  The  use  of  these 
constants  resulted  in  a  sputtering  ratio  versus  primary  energy  curve 
of  improper  shape;  the  curve  decreased  with  increasing  energy  with 
no  sign  of  the  experimental  Qj4J  maximum  around  3-  to  5-keV. 
Subsequently,  an  argon-copper  potential  was  employed  for  the  argon- 
copper  interactions,  while  the  Gibson  2  was  retained  for  the  copper- 
copper  system.  This  A-Cu  potential  is  also  of  the  Born-Mayer  form 
and  uses  semi-empirical  parameters  calculated  by  Harrison  ett  al. 
[26j  from  experimental  data  collected  by  Magnuson  and  Carlston  J33J  • 
The  sputtering  ratio  curve  obtained  with  this  two  potential  approach 
exhibited  a  peak  in  the  proper  region,  but  the  sputtering  ratio 
magnitude  was  too  low  at  high  energies.  To  increase  this  ratio  the 
A-Cu  potential  was  "hardened"  by  increasing  "A"  in  the  V(r)  expression 
from  the  calculated  value  of  11.435  to  12.25.  Figures  5,6  and  7 
result  when  the  tabulated  parameters  are  used  in  the  Born-Mayer 
potential  functions 


System 

A 

B  (lattice  unit) 

A-Cu 

12.25 

-7.597 

Cu-Cu 

10.024 

-9.197 

The  computer  program  is  a  mixture  of  Fortran  60  and  Fortran 
Symbolic  language  (see  appendix  II).  The  runs  were  conducted  on  a 
Control  Data  Corporation  1604  computer.  Running  time  for  each  impact 


point  was  five  minutes  for  the  (110)  sample,  nine  minutes  for  the 
(111)  and  from  ten  to  15  minutes  for  the  (100).  Very  generally,  the 
program  constructs  the  lattice,  bombards  the  crystallite  with  a 
single  primary,  performs  the  iteration  procedure  to  determine  forces 
and  displacements  at  the  proper  time  steps,  and  stops  when  the 
total  potential  energy  of  all  the  atoms  is  less  than  0.5  eV,  or  at 
a  time  step  chosen  by  experience  so  that  events  of  interest  have 
already  occurred.  The  output  consists  of  the  position  and  velocity 
of  each  atom  with  a  total  energy  greater  than  thermal,  assumed  to 
be  0.25  eV  (see  appendix  IV).  Output  is  received  when  the  run  is 
25$,  50$,  75%  and  100%  complete. 

At  the  end  of  the  runs  the  totals  of  kinetic  plus  potential 
energies  differed  from  the  primary  initial  energies  (the  lattice 
has  no  energy  until  disturbed)  by  an  average  of  1%,     This  error  is 
considered  well  within  the  uncertainties  inherent  in  the  model, 
especially  the  kinematic  simplifications  of  the  iteration  procedure. 


3.  Results 

A.  GENERAL  DESCRIPTION  OF  THE  SPUTTERING  PROCESS 
When  energy  from  the  primary  is  transferred  to  lattice  attorns 
through  the  n-body  collision  process,  these  atoms  (secondaries)  are 
displaced  from  their  sites.  The  secondaries  also  interact  in  a 
series  of  collisions,  and  eventually  some  will  gain  momentum  back 
towards  the  crystal  surface.  If  any  atoms  have  energies  great 
enough  to  overcome  the  surface  binding  energy  they  will  emerge  from 
the  crystal.  Very  generally,  this  is  the  sputtering  process. 

Cheney  etfc  al.  [jLJ  have  stated  that  sputtering  is  expected  to 
be  great  when  there  is  a  large  transfer  of  energy  from  the  primary 
to  the  crystal  lattice  occurring  reasonably  near  the  surface,  a 
rapid  reversal  of  the  momentum  vector  of  the  primary,  a  mechanism 
for  the  transfer  of  energy  back  to  the  surface,  an  easy  escape  of 
surface  atoms  and  a  comparatively  large  atomic  volume  or  mass. 
Relatively  open  crystal  structures  permitting  deep  penetration  of 
the  primary  results  in  a  low  sputtering  yield.  Thus,  the  yield  is 
dependent  upon  the  orientation  of  the  crystal  [24] •  However,  a 
crystal  which  has  a  history  of  radiation  damage  may  not  permit 
channeling  [20,32 1  in  open  crystallographic  directions,  in  which 
case  the  sputtering  yield  may  actually  increase.  Other  works  (_25,30J 
have  suggested  that  the  ejection  of  atoms  along  preferred  directions 
is  the  result  of  momentum  focusing  chains.  In  light  of  this 
research  some  of  these  concepts  must  be  modified. 


B.  THE  INDIVIDUAL  ORIENTATIONS 
1.  The  (110)  Sample 

Inspection  of  figure  8  reveals  that  most  of  the  ejected  atoms 
come  from  the  region  near  the  impact  rectangle.  A  detailed  analysis 
of  the  mechanism  whereby  a  few  of  these  atoms  are  sputtered  will 
lead  to  an  understanding  of  the  sputtering  process.  In  the  following 
discussion  "low"  energy  refers  to  1  keV  and  "high"  energy  refers  to 
10  keV  primary  energy. 

Atom  14  is  ejected  most  often.  At  low  energy  the  impact  points 
which  cause  sputtering  of  this  atom  are  those  in  the  upper  right 
quarter  of  the  representative  rectangle.  At  high  energy  the  impact 
points  in  the  center  of  the  rectangle  cause  it  to  sputter.  The 
number  of  impact  points  which  result  in  ejection  is  essentially 
independent  of  the  primary  energy. 

For  all  cases  in  which  this  atom  is  sputtered,  it  is  driven 
towards  the  crystal  interior  by  interaction  with  the  primary.  "Down" 
and  "left"  motion  is  also  imparted.  The  atom  is  repulsed  by  atoms 
lying  in  the  second  plane,  no.  33  usually,  and  is  reflected  towards 
the  front  surface.  This  process  imparts  momentum  to  no.  33  towards 
the  back  of  the  crystal,  which  accounts  for  the  relativly  few  times 
that  it  is  sputtered.  At  10  keV  no.  14  gains  enough  momentum  to 
interact  also  with  atom  no.  8  after  rebounding  from  no.  33. 

Atom  15  is  a  major  contributor  to  the  3  keV  peak  of  the  sputter- 
ing ratio  curve.  Impact  points  that  cause  this  atom  to  sputter  lie 
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in  the  right  center  portion  of  the  representative  rectangle  at  low 
energy,  and  the  bottom  center  at  high  energy.  No.  15  is  ejected 
less  often  with  increasing  primary  energy.  Channeling,  resulting 
in  deeper  ion  penetration  before  interaction,  occurs  with  greater 
primary  energy,  which  explains  this  behavior.  The  mechanism  which 
ejects  15  is  dependent  upon  interaction  at  the  front  surface.  The 
atom  is  initially  directed  to  the  right  by  the  primary.  The  ensuing 
collision  with  no.  16  rebounds  no.  15  out  of  the  crystal.  Addition- 
ally, at  higher  energies,  the  primary  rebounds  from  no.  14  and 
ejects  15  directly  by  getting  behind  it  before  interacting  strongly. 
Primaries  entering  at  impact  points  in  the  left  half  of  the  rectangle 
do  not  interact  strongly  enough  with  no.  15  to  impart  the  velocities 
necessary  to  collide  with  no.  16. 

Atom  37  is  also  a  contributor  to  the  sputtering  ratio  peak. 
The  impact  points  responsible  are  in  the  upper  left  quarter  of  the 
rectangle,  where  the  primary  "bypasses"  no.  14  and  interacts  mainly 
with  second  layer  atoms,  nos.  38  and  39  in  particular.  No.  38  is 
driven  in  and  to  the  left  in  a  direction  to  be  behind  no.  37  before 
approaching  closely  enough  to  ineract.  Thus  37  acquires  motion  out 
of  the  crystal  and  towards  the  left. 

Comparable  to  atom  no.  15,  the  impact  points  that  sputter  no. 
37  lie  in  a  region  where  the  primary  is  channeled  at  higher  energies, 
with  resultant  less  ion-atom  interaction  at  the  crystal  front  planes. 
No.  37,  then,  is  also  sputtered  less  at  higher  ion  energies. 
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At  low  energy  atom  19  is  initially  given  a  velocity  up,  in  and 
towards  the  left.  Interactions  with  no.  24  at  the  front  surface 
and  no.  43  in  the  second  plane  result  in  reflection  towards  the  front. 
The  impact  region  responsible  for  sputtering  no.  19  is  also  the 
upper  left  quarter,  with  a  few  points  around  the  perimeter  of  this 
area.  With  increased  primary  energy  impact  points  where  channeling 
occurs  cease  to  cause  sputtering  of  this  atom,  but  the  points  out- 
side of  this  area  still  eject  no.  19.  At  high  energy  the  primary 
is  reflected  from  the  third  layer  atoms  59  and  64,  and  takes  no.  19 
with  it  as  it  comes  back  out  the  front. 

At  10  keV  atoms  13,  38  and  40,  in  addition  to  nos.  14  and  15, 
are  sputtered  often.  Atom  13  is  sputtered  as  a  result  of  the  primary 
driving  14  behind  13,  whence  the  latter  acquires  motion  out  and  to 
the  left,  rebounding  from  no.  12  before  emerging.  Atom  38  is  sputter- 
ed by  an  identical  mechanism,  except  that  the  primary,  after  inter- 
acting with  no.  39  gets  behind  no.  38,  and  the  latter  interacts  with 
37  before  it  comes  out  of  the  crystal.  In  an  analogous  manner  no. 
39  sputters  atom  40. 

An  interesting,  but  not  surprising,  phenomenon  is  the  number 
of  atoms  sputtered  from  deep  (third  or  greater)  layers  as  a  function 
of  primary  energy.  There  are  none  sputtered  at  1  keV,  one  at  3  keV, 
four  at  5  keV,  and  six  at  10  keV.  This  behavior  would  be  expected 
since  the  primary  penetrates  further  with  increased  energy  and 
transfers  momentum  to  the  deeper  atoms. 
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As  mentioned  previously  in  section  2,  the  two  potential  approach 
resulted  in  a  sputtering  ratio  curve  of  the  proper  shape.  To  adjust 
this  curve  to  the  ordinates  obtained  experimentally  Q2^H ,  a  "thresh- 
old energy"  notion  is  employed.  An  arbitrary  value  is  chosen,  and 
the  assumption  made  that  any  atom  with  kinetic  energy  less  than  this 
amount  (in  the  direction  normal  to  the  surface)  will  not  escape  the 
crystal.  A  family  of  curves  result  (see  fig.  5)  from  which  an 
estimate  can  be  made  of  the  binding  energy  of  atoms  to  a  crystal 
surface.  Thus  4  eV  is  the  value  estimated  for  the  copper  (110) 
surface  since  this  figure  yields  a  curve  with  the  best  fit  to 
experimental  data. 

An  analysis  of  the  angular  distribution  of  the  sputtered  atoms 
leads  to  some  extremely  interesting  results.  At  1  keV  primary 
energy  one  third  of  the  ejected  atoms  leave  the  surface  with  equal 
velocities  in  the  "x"  and  nzw  directions.  These  velocities  range 
from  one-half  to  equal  the  "y"  velocity.  That  is,  a  large  percent- 
age of  atoms  emerge  approximately  in  the  ^111^  direction,  but  nearer 
to  the  surface  normal.  One-sixth  of  the  atoms  emerge  in  the  ^110^ 
direction  while  the  remainder  are  randomly  oriented.  Thus  the  model 
predicts  a  deposit  spot  pattern  consisting  of  a  central  spot,  four 
^111^  spots  and  considerable  haze. 

At  10  keV  one-third  of  the  atoms  emerge  in  the  ^110^  direction, 
one-fifth  in  approximately  the  ^111^  direction,  and  one-fifth  near 
the  ^100^  direction.  At  this  energy  the  model  predicts  the  central 
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spot,  four  ^111^  spots,  plus  two  ^100^  spots  which  may  merge  with 
the  central  spot,  and  less  haze  than  obtained  at  1  keV. 

2.  The  (111)  Sample 

As  an  example  of  the  sputtering  mechanism,  consider  atom  19, 
which  is  ejected  most  often.  This  atom  has  motion  imparted  by  the 
primary  "up'1  and  towards  the  right  (see  fig.  3).  The  collision 
with  atom  25  alters  its  momentum  vector  towards  the  front  surface. 

The  sputtering  ratio  curve  is  shown  in  figure  7.  Although  the 
shape  of  the  curve  approaches  the  experimental  data  for  6  eV  threshold 
energy,  the  magnitudes  disagree  significantly.  Inspection  of  the 
figure  shows  that  different  threshold  energies  do  not  yield  a  family 
of  curves,  but  that  the  curve  shape  is  altered  with  changes  in 
energy. 

The  deposit  pattern  for  this  orientation  consists  of  a  central 
spot  with  three  large  surrounding  spots  corresponding  to  the  ^011^ 
directions. 

3.  The  (100)  Sample 

Analysis  of  the  paths  of  sputtered  atoms  reveals  that  the 
ejection  mechanism  is  consistent  with  the  other  two  orientations. 
Figure  9  shows  that  the  atoms  near  the  impact  area  are  ejected  most 
often. 

The  sputtering  ratio  curve  is  shown  in  figure  6.  There  is  close 
agreement  with  the  experimental  data  L.24J  ^   tne  1~5  keV  range.  At 
10  keV  the  sputtering  ratio  is  too  high,  indicating  that  the  potential 
function  may  be  too  hard  at  higher  energies. 
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The  deposit  pattern  for  this  orientation  consists  of  a  heavy 
concentration  of  atoms  in  the  center,  with  lesser  concentrations 
in  diagonally  opposite  quadrants. 
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4.  Discussion 

A.  The  (110)  Surface 
As  mentioned  in  section  3,  the  majority  of  the  atoms  which  are 
sputtered  come  from  the  region  around  the  impact  area.  This  obser- 
vation indicates  a  model  which  accounts  for  sputtering  as  the  result 
of  the  direct  transfer  of  momentum  from  the  primary  ion  to  the  lattice 
atoms,  rather  than  by  chaining.  The  mechanism  whereby  crystal  atoms 
are  ejected  supports  this  view.  Although  there  are  differences  in 
this  mechanism  at  high  and  low  energies,  there  is  no  evidence  in 
this  study  that  momentum  focusing  plays  a  role  in  the  sputtering 
process. 

The  threshold  energy  necessary  to  achieve  agreement  with  the 
experimental  sputtering  ratio  curve  is  considered  to  be  quite  reason- 
able. The  value  of  4  eV  is  typical  of  the  heat  of  sublimation,  which 
is  less  than  the  energy  required  to  displace  an  atom  from  its  lattice 
site.  The  displacement  energies  which  have  been  determined  are 
applicable  to  interior  atoms,  while  surface  atoms  are  expected  to 
be  less  tightly  bound. 

The  model  successfully  predicts  the  shape  of  the  sputtering  ratio 
curve.  Discrepancies  between  experimental  and  theoretical  magnitudes 
can  be  attributed  to  the  approximations  of  the  potential  function 
in  the  theory  or  the  modification  of  the  crystal  during  irradiation 
experiments • 

The  deposit  spot  pattern  observed  by  Nelson,  Thompson,  Farmery 
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and  Hall  [37J  is  also  successfully  predicted.  Although  the  experiment- 
al group  was  working  at  higher  energies,  deposit  patterns  above  10 
keV  are  expected  to  have  the  same  general  features.  At  1  keV  this 
study  predicts  a  spot  pattern  similar  to  those  obtained  by  Southern 
et^  al.  ("25J.  This  group  worked  with  a  ^01]^  sample  but  crystall- 
ography considerations  show  that  the  ^110^  should  have  similar 
patterns.  The  model  agrees  qualitatively  with  the  experimental 
results  of  the  latter  group  in  two  important  respects;  (1)  the  non- 
central  ejection  spots  are  closer  to  the  surface  normal  than  would 
be  expected  by  consideration  of  crystal  geometry,  and  (2)  the  neigh- 
boring rows  influence  the  motion  of  the  ejected  atoms. 

B.  The  (111)  Surface 

The  sputtering  mechanism  is  similar  to  the  (110)  sample  in  that 
momentum  focusing  is  not  evident.  The  atoms  are  ejected  as  the 
result  of  sequences  of  first  and  second  layer  collisions. 

The  estimated  binding  energy  is  taken  to  be  6  eV,  the  value 
which  yields  the  proper  shape  of  the  sputtering  ratio  curve.  Dis- 
crepancies in  magnitude  between  the  model  and  experiment  are  consider- 
ed to  be  the  result  of  failure  of  the  potential  function  parameters. 
"Hardening"  the  potential  for  future  studies  of  this  orientation  may 
give  better  agreement. 

Selection  of  different  binding  energies  have  the  greatest  effect 
on  the  magnitude  of  the  sputtering  ratio  at  lower  energies.  This 
indicates  that  lower  primary  energies  result  in  ejection  of  lower 
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energy  atoms. 

The  deposit  pattern  experimentally  observed  by  Southern  et„  al. 
[25j  is  successfully  predicted. 

C.  The  (100)  Surface 

Results  of  bombardment  of  this  surface  also  indicate  that 
momentum  focusing  is  not  evident.  Sputtering  is  the  result  of  first 
or  second  layer  collisions. 

The  estimated  binding  energy  is  2.6  eV.  Additional  investigation 
in  the  5-10  keV  range  is  being  conducted  to  bring  the  sputtering 
ratio  curve  into  closer  agreement  with  experimental  observations. 
The  shape  of  the  curve  is  sensitive  to  the  threshold  energy,  which 
has  the  greatest  effect  on  magnitude  at  3  keV  (see  fig.  6).  The 
possibility  exists,  in  nature,  that  the  binding  energy  of  the  crystal 
is  a  function  of  the  amount  of  disturbance,  and  hence  a  function  of 
primary  energy. 

The  deposit  pattern  corresponds  to  "half"  of  the  pattern 
observed  by  Southern  et  al.  |~25 J.  That  is,  there  is  a  central  spot 
and  two  ^011>  spots.  Bombarding  the  crystallite  in  a  mirror  image 
impact  triangle  may  yield  the  other  two  ^Olly  spots. 


18 


5.  Conclusions 

This  computer  simulation  program  is  considered  adequate  in  most 
respects.  The  model  predicts  the  proper  shape  of  the  sputtering 
ratio  curves,  the  angular  distribution  of  ejected  atoms,  binding 
energies  in  the  neighborhood  of  expected  values,  and  magnitudes  of 
the  sputtering  ratio  for  two  of  the  three  orientations. 

Momentum  focusing  does  not  enter  into  the  sputtering  mechanism 
for  the  argon-copper  system  in  the  1-  to  10-keV  range.  However, 
there  are  cases  of  chain  transmission  of  momentum  of  surface  atoms 
parallel  to  the  surface. 

There  is  a  wealth  of  experimental  data  available  which  requires 
theoretical  exploration.  Further  work  using  this  model  should 
include  study  of  the  energy  distribution  of  ejected  atoms,  sputtering 
ratio  as  a  function  of  primary  incoming  direction,  inclusion  of 
thermalized  motion  of  lattice  atoms,  and  investigation  of  transmission 
sputtering,  all  for  additional  systems  of  elements  and  other  energy 
ranges.  The  computer  program  should  be  expanded  to  include  a  graph 
plot  output  of  the  sputtering  ratio  curve  and  a  plot  of  the  angular 
and  energy  distributions  of  the  ejected  atoms. 

The  computer  programming  assistance  of  LT  Peter  D.Stogisis 
gratefully  acknowledged.  I  am  indebted  to  Prof.  Don  E.  Harrison  Jr. 
for  his  guidance  and  assistance  in  all  phases  of  the  work. 

Note  in  support  of  argument.  Recent  work  on  argon  sputtered 
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lithium  by  Stein,  Malakhof  and  Smith  [38J  led  the  authors  to 
suggest  that,  at  low  energy,  there  is  small  likelihood  of  extensive 
penetration  of  the  lattice  by  the  primary.  Hence,  there  is  little 
chance  for  a  large  number  of  n-body  collisions  yielding  a  focused 
momentum  change  to  take  place.  The  lobes  of  the  distribution 
pattern  of  the  sputtered  atoms  are  believed  to  be  created  by  a 
limited  number  of  n-body  collisions  near  the  surface.   (This 
effect  should  be  even  more  pronounced  for  copper  since  the  mass 
of  the  target  atoms  is  considerably  greater  than  lithium) . 
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FIG.   1 
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FACE  CENTERED  CUBIC  (100) 
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FIG.  2 
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Face  Centered  Cubic    (III) 


FIG.   3 
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FACE  CENTERED  CUBIC  (110) 


FIG. 8 
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APPENDIX  I 


CRYSTAL  GEOMETRY 


A,  Lattice  Dimmensions 


Orientation 

Dimensions  (L.U.) 

X 

Y 

Z 

(no) 

9  x  .707  =  6.36 

4  x  .707  =  2.83 

8  x  1.0  =  8.0 

(100) 

6  x  1.0  -  6.0 

4  x  1.0  =  4.0 

8  x  1.0  =  8.0 

(111) 

9  x  .707  =  6.36 

4x1.155  =4.62 

5.67x1.22=  6.94 

1  L.U.  =  1.807 
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B.  Atomic  Numbering  System 


Y-plane  number 

Orientation 

1 

2 

3 

4 

5 

(no) 

2-26 

27-46 

47-71 

72-91 

92-116 

128-159 

122-151 

(100) 

2-33 

34-64 
32-61 

65-96 

97-127 

(111) 

2-31 

62-91 

91-121 
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APPENDIX  III 

PROGRAM  DISCUSSION 

A.  Brief  Explanation  of  the  Program. 

The  descriptions  below  pertain  to  the  sections  of  the  program  as 

outlined  in  appendix  II. 

Section  Is  Allocates  space  in  computer  cells  for  arrays  of  various 
properties  of  the  atoms;  establishes  input  and  output 
formats;  establishes  the  Born-Mayer  form  for  both 
potentials. 

Section  2  s  Defines  various  constants  and  variables;  reads  input  data. 

Section  3s  Constructs  the  crystallite. 

Section  4s  Commences  the  bombardment  by  assigning  the  primary  its 
proper  direction  and  velocity. 

Section  5s  Force  double  iteration  procedure.  This  section  calculates 
the  forces  on  the  atoms  by  use  of  a  triangular  matrix, 
then  later  recalculates  the  forces  after  the  atoms  have 
been  displaced  by  the  original  forces.  The  average  of 
these  forces  is  then  applied  to  the  initial  positions 
of  the  atoms  to  determine  their  displacements. 

Section  6s  Determines  the  tentative  positions  of  the  atoms  under 
the  forces  calculated  in  section  5,  then  returns  the 
program  to  section  5  for  force  recalculation. 


44 


Section  7s  Establishes  positions  of  the  atoms  under  the  average  force. 
Calculates  velocities  of  all  atoms. 

Section  8s  Tests  the  current  time  step  to  see  if  output  is  called 
for.  If  affirmative,  prints  output  headings. 

Section  9s  Calculates  atomic  separations  from  the  positions  deter- 
mined in  section  7.  Calculates  potential  energies  based 
upon  these  separations  and  the  Born-Mayer  parameters. 

Section  10s  Computes  kinetic  energies  of  the  atoms  from  their  velocities. 
Calculates  the  total  binding  energy  of  escaping  atoms. 

Section  lis  If  the  time  step  is  such  that  output  is  called  for  this 
section  prints  displacement,  velocity,  kinetic  and 
potential  energies  of  each  atom  with  a  total  energy 
greater  than  0.25  eV. 

Section  12s  Indexes  any  atoms  which  have  escaped  from  the  crystal  so 
that  they  will  be  excluded  from  future  "do"  loops  (an 
attempt  to  reduce  computer  running  time).  Recalculates 
the  time  step  interval  on  the  basis  of  the  most  energetic 
atom  still  in  the  crystal.  Tests  the  total  potential 
energy  and  time  step  interval  against  values  selected 
to  be  indicative  of  the  end  of  a  run.  Continues  the  run 
if  tests  are  negative. 

Section  13s  Prints  the  status  report  at  the  end  of  the  run,  which 
includes  the  position,  velocity,  kinetic  and  potential 
energies  of  each  atom.  Returns  the  program  to  the  input 
step  for  reception  of  the  next  impact  point  data. 
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B.  Terms  not  defined  in  Leeds*  Program  (_34J. 

PPOTF(X) .PFORCF(X) ,PFXA,PFRC  and  PPTC  are  the  argon-copper  functions 

which  correspond  to  the  copper-copper  functions  POTF(X),FORCF(X), 

FXA,FRC  and  PTC. 

LCUT  indexes  any  atoms  which  have  been  ejected  from  the  crystal  as 

explained  under  section  12,  appendix  II  A. 

NLDTI  is  the  number  of  time  steps  between  recalculation  of  DT,  the 

time  step  interval. 

NL  is  an  index  used  to  test  if  NLDTI  time  steps  have  elapsed  since 

DT  was  last  recalculated. 

QUIT  is  the  value  of  total  potential  energy  chosen  to  test  for 

program  shutdown. 

SLOW  is  the  value  of  the  time  step  interval  chosen  to  test  for 

program  shutdown. 

BENERGY  is  the  binding  force  assigned  to  the  crystal  surface,  in  eV. 

CVB  is  the  conversion  factor  to  make  BENERGY  consistent  with  units 

employed  within  the  program. 

BRY  is  the  cumulative  distance  from  the  front  surface  of  all  atoms 

that  have  been  ejected. 

BF  is  the  cumulative  binding  energy  of  these  ejected  atoms.  The  value 

is  determined  by  a  force  x  distance  calculation  (BENERGY  x  CVB  x  BRY). 
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